Nanotubular Boron-Carbon Heterojunctions 
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Linear nanotubular boron-carbon heterojunctions are systematically constructed and studied with 
the help of ab initio total energy calculations. The structural compatibility of the two classes of ma- 
terials is shown, and a simple recipe that determines all types of stable linear junctions is illustrated 
in some detail. Our results also suggest the compatibility of various technologically interesting types 
of nanotubular materials, leading to novel types of nanotubular compound materials, and pointing 
out the possibility of wiring nanotubular devices within heterogeneous nanotubular networks. 

PACS numbers: 81.07.Lk, 81.07.Dc, 73.22.-f, 61.46.+W 



INTRODUCTION 

Carbon nanotubes (CNTs) Q have a number of fasci- 
nating properties, which could make them the material 
of choice for the future miniaturization of electronic de- 
vices: They are good thermal conductors, exhibit a good 
thermal resilience, are ductile and superb field emitters 
Q. Furthermore their electronic properties range from 
metallic to semi-conducting, depending on their chirali- 
ties and radii |2j- 

Besides there are many other nanotubes made from 
inorganic materials 0, , like BN [(| or M0S2 nan- 
otubes which have already been synthesized, or materials 
that have been predicted by theory like metal-boron nan- 
otubes H 0| or the structurally related CaSi2 nanotubes 
|lf)j ] . Another promising new type of nanotubular materi- 
als are nanotubes consisting of pure boron (B-nanotubes) 
[ill fl2| . The stability and the mechanical properties of 
B-nanotubes (BNTs) would be quite similar to C- and 
BN-nanotubes, but all BNTs should be metallic, inde- 
pendent of their chirality Very recently Ciuparu 
et al. 14] successfully synthesized these nanotubes and 
thus confirmed the existence of BNTs, after similar ef- 
forts had already lead to the discovery of novel types of 
boron nanowires by various other groups 0, 0, LL3, 0] . 

With such a large number of nanotubular systems 
available (a number which is very likely to further in- 
crease in the near future), possible nanotechnologies are 
most likely to be based on this variety of nanotubular ma- 
terials rather than just CNTs (see [ljj,l20]). The function- 
alities of such heterogeneous nanotubular networks would 
mainly arise from their composite character. 

But before one can seriously start any discussion about 
the future impact of heterogeneous nanotubular net- 
works, it will be necessary to explore the stability of 
their key elements, first: the interfaces between differ- 
ent nanotubular materials. In the following, we will, for 
the first time, study the compatibility of different nan- 
otubular materials in a systematic fashion by examining 
the structural and electronic properties of suitable model 
junctions between those systems. 

To this end we decided to focus our studies on boron- 



carbon nanojunctions, which may be seen as a struc- 
tural paradigm for a large class of junctions, compris- 
ing most of the non-carbon based nanotubular materials 
mentioned above. On the basis of three exemplary struc- 
tures, we will demonstrate the compatibility of CNTs 
and BNTs. Due to the exemplary character of the model 
junctions chosen for this study, our results also suggest 
the existence of similar junctions between a large number 
of related nanotubular materials. 



NANOTUBULAR HETEROJUNCTIONS 

Nanotubes are geometrically constructed by rolling up 
a rectangular sheet that has been cut from a planar struc- 
ture. For CNTs this planar structure will be the hon- 
eycomb lattice [jj|, while for BNTs the planar reference 
structure will be a (puckered) hexagonal lattice [Till ]. 

In order to approach the rather complex problem of 
heterogeneous nanotubular junctions we restricted our 
examinations to linear B-C junctions with constant chi- 
rality, where both parts of the junctions have similar 
radii. The sought B-C nanojunctions (and all other 
structures considered in this paper) are directly related 
to the structure of pure CNTs, as they can obviously 
be constructed from a single sheet. Therefore we may 
also classify them in a standard fashion employing a 
pair of integers (N,M) which specify the chiral vector 
Ch = Nb\ + Mb-2 . Here bi and b% are basis vectors for 
a honeycomb lattice (see Fig^). 

Now the general prescription to create our sample B- 
C junctions goes as follows: A rectangular sheet is cut 
from a honeycomb lattice, such that its horizontal sides 
equal the chiral vector, and the vertical sides (dubbed 
^-direction in the following) being even multiples of a 
basic period, depending on the chirality of the tube 0]. 
An imaginary horizontal line in the middle of the sheet 
separates the upper and the lower part. The lower part 
consists of the original honeycomb lattice and is deco- 
rated with carbon atoms at each lattice site. The upper 
part is transformed into a hexagonal lattice after insert- 
ing another lattice point at the center of the honeycombs. 
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FIG. 1: Basic construction of B-C junctions, (a): (9,6) sheet, units are A, light atoms are boron, dark atoms are carbon. Also 
indicated are the chiral vector Ch and the basis vectors of a honeycomb lattice 61,62- The encircled parts emphasize the two 
structure elements A and B that appear in the transition region, and the dashed line separates the boron part from the carbon 
part, (b): chiral nanojunction constructed from (a) by rolling up the sheet and gluing it at the vertical sides, where the grey 
cylinder should emphasize the resulting tubular structure. 



Next, this structure is decorated with boron atoms. Fi- 
nally the sheet is rolled up and glued along its vertical 
sides (see Fig. P). 

For the stability of these junctions it will be crucial to 
understand the local atomic configurations close to the 
interface. Using the generation scheme described above, 
only two different structure elements, dubbed A and B 
in the following, may form close to the transition region 
(emphasized in Fig. QJ. Thus a detailed study of the 
basic structural and chemical properties related to the 
formation of A and B will be sufficient to understand the 
stability of any linear and chiral nanojunction. These 
structure elements located around the interface regions 
will imply a planar coordination number of four for the 
constituent carbon atoms. Such coordinations are not 
found in pure CNTs, but higher planar coordinations of 
carbon atoms in contact with boron atoms have already 
been reported by Exner et al. |2l| for small BC-clusters. 
We also note that non-chiral armchair and zigzag systems 
contain only one type of structure element A or B, and 
in order to be able to study their effects in a systematic 
fashion, we will restrict ourselves to armchair and zigzag 
types of nanojunctions. 



THEORETICAL METHODS 

Due to its electron deficient character [2^] boron has a 
complicated and versatile chemistry. The only theoretical 
tools that allow to describe its chemistry properly are 
first principles calculations |23|. 

In order to carry out ah initio type of simulations for 
nanotubular compound systems, we have to construct 
a solid composed of suitable unit cells containing those 
junctions. Within the xy-plane, we arrange the tubular 
junctions side by side on a hexagonal lattice of lattice 
constant a. Therefore we will simulate bundles (ropes) 
of linear junctions rather than single free-standing nano- 
junctions. In the z-direction we simply pile up our unit 
cells (with lattice constant c). The whole procedure will 
lead to tubular structures that are a linear array of B-C 
and C-B transitions (see Fig. 0). Due to the periodic 
boundary conditions these systems can at best approxi- 
mate single nanojunctions that would be present in the 
limit a, c — > 00. Nevertheless the small approximants 
considered in this paper are sufficient to understand the 
properties of the structure elements A and B, as shown 
below. 

The numerical calculations themselves were carried out 
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FIG. 2: The unit cell for hexagonal tubular bundles (ropes), 
with its lattice constants a and c, and the atomic decoration 
for the unrelaxed approximant Zl (see text). Light atoms are 
boron, dark atoms are carbon and the grey cylinder again 
emphasizes the tubular structure of the junction. 

using the VASP ab initio package, version 4.4.5 [24II25I . 
The latter is a density functional theory (DFT) |26J 
based ab initio code using plane wave basis sets and a 
(super) cell approach to model solid materials. During 
all simulations, the electronic correlations were treated 
within the local-density approximation (LDA 27]), and 
the ionic cores of the system were represented by ultra- 
soft pseudopotentials |28|. With the help of the VASP 
program, one can determine interatomic forces and relax 
all degrees of freedom for a given decorated unit cell in 
order to detect atomic configurations which correspond 
to (local) minima on the total energy hypersurfaces. In 
order to perform these kinds of structure optimizations 
we employed a preconditioned conjugate gradient algo- 
rithm [29j and let all degrees of freedom relax, i.e. the 
complete set of atomic configurations as well as the unit 
cell parameters. The total energy as well as the fc-point 
sampling were converged such that changes in the total 
energies were less than 10 -3 eV and interatomic forces 
were less than 0.05 eV/A. 

It turned out that the energy hypersurfaces of our 
structures were rather complex and full of local minima. 
In order to approach the most stable structures, we devel- 
oped a procedure where the starting configurations were 
prerelaxed with lower numerical precision until a (local) 
minimum was found, and afterwards we continued the 
relaxations with optimal precision. A reduced precision 
leads to somewhat imprecise interatomic forces, but we 
found that such a procedure would result in the scanning 
of the energy hypersurface over a wider range. 

In the following, we will discuss our results for the 
following types of model junctions: Armchair 1 (Al) and 
Armchair 2 (A2) approximate an armchair (6,6) junction, 
and Zigzag 1 (Zl) approximates a zigzag (6,0) junction 
(Fig. EJ). Besides (6,6) and (6,0) BNTs and CNTs have 
been relaxed, which are considered as suitable reference 
systems (Fig. |3Jl. 
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(a) Armchair (6,6) CNT (b) Armchair (6,6) BNT 
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(c) Zigzag (6,0) CNT (d) Zigzag (6,0) BNT 
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FIG. 3: Top view of the (relaxed) reference structures placed 
around the center of their unit cells. Encircled atoms in (b) 
have additional covalent bonds to neighboring nanotubes on 
the superlattice. Units are A, light atoms are boron and dark 
atoms are carbon. 



Al was prerelaxed with the Brillouin zone being sam- 
pled by just the central T-Point, and finished with a 
4x4x4 grid. A 4x4x4 grid was used to finish Zl, with 
smaller FFT-grids and a smaller plane wave basis set to 
do the prerelaxation. Finally A2 was prerelaxed with the 
T-Point, only and finished with a 2x2x2 grid. The (6,0) 
CNT was relaxed with a 4x4x7 grid and the total en- 
ergy was calculated on a 5x5x7 grid. The armchair (6,6) 
CNT was prerelaxed with 4x4x6 and finished with 5x5x9. 
The boron (6,0) zigzag nanotube was started with 3x3x5 
and finished with 5x5x11 fc-points. The relaxed arm- 
chair BNT structure was obtained from || , where it was 
relaxed on a 7x7x7 grid. The cutoff energy for the ex- 
pansion of the wave function in plane waves was 286.6 
eV for A2, 358.2 eV for Al, Zl and the CNTs, and 321.4 
cV for the BNTs. 

The cohesive energy of each structure (Tab. was 
calculated by dividing the binding energy per unit cell 
by the number of atoms contained within that unit 
cell. In order to judge the stability of our model junc- 
tions, we decided to compare the cohesive energies of 
the relaxed structures to a reference cohesive energy 
of a phase-separated system of pure CNTs and pure 
BNTs of the same structure type, see Table [I] For 
the (6,6)-junction the reference binding energy would be 
E^ d = 24 x 10.023 + 36 x 7.004 = 492.696 eV, and 
the corresponding reference cohesive energy would be 
E r c l{ = E T h f nd / 60 = 8.212 eV/atom. All cohesive ener- 
gies as well as energetic differences with respect to the 
reference cohesive energies are given in table H| 
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TABLE I: Structure data and energies for nanojunctions and reference structures. (N,M): nanotube structure type, C: number 
of carbon atoms per unit cell, B: number of boron atoms per unit cell, (a,c): lattice constants a and c of unit cell in A, ac-c, 
(ib-b, csb-c: range of bond lengths in A of C-C, B-B, and B-C bonds, respectively, E co ^\ cohesive energies in eV/atom, E^: 
reference cohesive energy in eV/atom (see text). 
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. 1.78 
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"relaxed structure from 0] 



RESULTS 

Our survey started with the relaxation of the refer- 
ence structures, which are bundles (ropes) of pure CNTs 
and BNTs. They are displayed in Fig. [3] (top view). 
While the surface of small radius CNTs is always smooth 
(Fig. Ek,c), BNTs can exhibit a puckered surface [l3j . 
which shows up in Fig. |3t>, but not in Fig. Due to 
higher B-B bond lengths, the average radius of a BNT 
is always bigger than the radius of a CNT of the same 
structure type. Another major difference between CNTs 
and BNTs is the potential of the latter to form covalent 
intertubular bonds Q- In contrast, the CNTs may only 
bind to each other via van der Waals types of interac- 
tions Such covalent bonds are indeed present in the 
model (6,6) BNT, and the atoms that form these bonds 
are encircled in Fig. |3Jd. The coordination number of 
boron atoms in BNTs is usually six t 3Q| or even seven 
for atoms forming intertubular bonds, and it is always 
three in CNTs. The range of bond lengths, the parame- 
ters of the unit cells, and the cohesive energies of all fully 
relaxed structures can be found in table [I] 

The relaxed approximants Al, A2, and Zl as well as 
the initial configurations are shown in planar view in Fig. 
0J The left column shows the sheet, from which we con- 
structed the starting configurations. In the middle the 
optimized junctions were projected back onto a similar 
sheet. The right column shows the junctions and their 
unit cells in top view. All of these structures have a 
Cg-Symmetry, e.g. they may be dissected into six iden- 
tical 60°-wedges (indicated by dashed lines in the planar 
views) j3l| . 

The cohesive energies of all approximants are below the 
cohesive energies of the corresponding phase-separated 
reference systems (see Tab. HJ. But the systems of iso- 
lated BNTs and CNTs are only sightly more stable than 
the junctions presented in this study. Each junction on 
the superlattice will bind to neighboring tubes via co- 
valent boron atoms. There are no covalent intertubular 
bonds involving carbon. Instead, the carbon part of the 



junctions almost rigidly stays in its initial geometry, and 
generally leads to smooth surfaces for the carbon part of 
the junctions. The boron on the other hand turns out to 
be more flexible, always forming puckered surfaces. 

We studied the properties of structure element A by 
setting up the approximants Al. The cohesive energy of 
Al lies 0.13 eV/atom below its cohesive reference energy, 
making this system slightly less stable than the phase 
separated reference system. After relaxation the hexag- 
onal symmetry of the boron part has disappeared, and 
the junction was transformed into a new structure with 
4-fold intratubc coordination for boron, and 3-fold coor- 
dination for carbon (see Fig^J)). It seems that carbon 
tends to retain the planar sp 2 -environment and forces the 
boron atoms into an unusual quasiplanar arrangement 
with rather low coordination numbers. 

To overcome the somewhat constrained arrangements 
of Al, we set up another approximant (not shown) cho- 
sen to be twice as long in the z-direction, and compris- 
ing twice as many atoms, but still containing structure 
element A, only. The relaxed structure had the same 
structural problems as Al and was almost isoenergetic to 
Al. This gives a consistent picture and certainly shows 
that the structural problems encountered for these types 
of armchair approximants are no artefacts of small unit 
cells. Thus the starting geometries involving the struc- 
ture element A are obviously far from optimal. 

In a next step we studied the structure element B by 
surveying approximant Zl. Its cohesive energy lies just 
0.07 eV/atom below the reference energy. The planar 
view in Fig. 0^ looks almost identical to its starting con- 
figuration and the basic lattice structures in both parts 
of the junction are maintained, allowing for high coordi- 
nation numbers of the constituent atoms. It seems as if 
this type of heterojunction could be formed by sticking 
together two perfect zigzag tubes. It therefore looks as if 
nanotubular junctions formed with structure element B 
are generally favorable. 

Therefore we tried to incorporate B in junctions of 
armchair type as well. To this end, one has to mod- 
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(a) Al - unrelaxed (b) Al - relaxed (c) Al - top view 




(d) Zl - unrelaxed (e) Zl - relaxed (f) Zl - top view 




(g) A2 - unrelaxed (h) A2 - relaxed (i) A2 - top view 




FIG. 4: First column: starting geometries in planar view, encircled parts emphasize the structure elements A and B. Dashed 
lines indicate six similar wedges partitioning the unit cells. Second column: relaxed structures in planar view. The local density 
of states of the encircled atoms in (e) and (h) is plotted in Fig. 03 Third column: top view of relaxed structures in their unit 
cells. Units are in A, and light atoms are boron, and dark atoms are carbon. 



ify the construction scheme and allow for a mixing of the 
species at the interface. Such a scenario has been realized 
in approximant A2 (Fig. third line) . The unit cell has 
been shifted by c/4. Now the two transition regions in 
every unit cell are clearly visible. The relaxed junctions 
shows properties very similar to Zl. Its cohesive energy 
is 0.07 eV/atom below the reference energy. Again the 
planar view in Fig. 0Ji looks almost identical to its initial 
configuration, indicating a starting point on the energy 
hypersurface very close to the sought minimum. Though 
very different in geometry Zl and A2 exhibit very similar 
properties: the same stability, the same mean coordina- 
tion number and the same mean number of intertubular 
bonds. 

Finally we would like to discuss the basic electronic 
properties of Zl and A2. The first line of Fig. [5] shows 
the total electronic density of states (DOS) of the struc- 
tures near the Fermi energy Ep. The double-plots al- 
low for a comparison of the local DOS of a central C- 
or B-atom (these atoms are encircled in Fig. 0^ and 
Eh)H2 to the DOS of the reference systems, which are 
considered as bulk. Both are very different. Thus the 



approximants obviously do not contain bulk states and 
the transition to the bulk proceeds over wider ranges, 
which cannot be simulated by the small approximants 
presented in this study. The position of major van-Hove- 
peaks and band gaps in the local DOS is homogeneous 
throughout the junction (thus it does not differ in the 
local and total DOS). This could be an indication for 
a common 7r-electron system extending throughout the 
junction. Finally we also found that the total DOS looks 
more similar to the local DOS of the boron atom rather 
than to the local DOS of carbon. 



CONCLUSIONS 

For the first time we explored the stability, geometry 
and basic chemical and electronic properties of sample 
carbon-boron heterojunctions. 

Starting from a general prescription to generate linear 
nanotubular junctions with arbitrary radius and chiral- 
ity we showed that only two structure elements A and 
B can appear at the B-C transition of such a junction. 
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FIG. 5: Total and local density of states (DOS) in arbitrary 
units of Zl and A2 (first and second column, respectively). 
In the split-panels we compare the total DOS of the reference 
structures (bulk) (see Fig. [3J to local DOS of atoms within 
the junctions. These atoms are encircled in Fig. 0^ and|lji. 
E — equals the Fermi energy Ef- 



The general atomic arrangements of A and B look rather 
similar, but in A the central carbon atoms will bind to 
two boron and two carbon atoms, whereas for B the cen- 
tral carbon atom will bind to three boron and one carbon 
atom, only. 

We then studied the properties of A and B indepen- 
dently by examining approximants for a (6,6) armchair 
and a (6,0) zigzag junction. From a structural point of 
view, the carbon part of all relaxed junctions looks rather 
similar to what is known for stand-alone carbon nan- 
otubes - exhibiting a smooth surface and preferring car- 
bon in a 3-fold coordinated environment [3j. The boron 
part of the junctions is always puckered and forms in- 
tcrtubular bonds, as known from stand-alone boron nan- 
otubes 0, |3(| , but the ideal hexagonal arrangement in 
the boron part of the junctions is only preserved by junc- 
tions that feature the structure element B. Such sys- 
tems have similar properties; they display the same aver- 
age coordination numbers, the same types of intertubu- 



lar bonds, and they are only 0.07 eV/atom less stable 
than a phase separated system of carbon and boron nan- 
otubes. Junctions that feature element A have structural 
problems and their formation energies per atom are ap- 
proximately twice as high as for junctions formed with 
element B. 

The central carbon atom of structure element B is 4- 
fold coordinated. Besides the known tetragonal sp 3 ar- 
rangements, we found new planar and quasiplanar atomic 
arrangements induced by the electron deficiency of boron 
|22| . Such forms have not been reported in solid system 
before and are in partial agreement with results found by 
Exner et al. [2l| , who studied small planar BC-clusters. 

The analysis of the total and local electronic DOS in- 
dicates that all junctions consist of interface states, only. 
A transition towards the bulk will proceed over wider 
ranges than those that were studied here. Furthermore 
there are indications for a common 7r-electron system ex- 
tending throughout the whole junction. 

After relaxation all approximants were still tubular, 
had distinct boron and carbon parts and showed a clear 
B-C interface. This itself is a strong indication for the 
compatibility of these two nanotubular materials. In de- 
tail, it seems that structure element B is the key to gen- 
erate optimal zigzag and armchair junctions, and those 
results also suggest that any stable linear and chiral nan- 
otubular B-C junction may be formed by building struc- 
ture elements B within the transition region. The relative 
instability of 0.07 eV/atom for the apparently optimal 
model junctions is surprisingly small, considering that 
at the interface both the boron and the carbon lattices 
are significantly disturbed, and second, that our approx- 
imants are really small. This energy barrier can be over- 
come thermodynamically, and it is very likely to decrease 
further for bigger systems. The influence of lattice per- 
turbations is apparent in Al, where only small disloca- 
tions of the boron atoms give rise to a doubling of the rel- 
ative instability. We thus conjecture that heterogeneous 
junctions between boron and carbon nanotubes should 
be possible. Our study also clarifies the local chemistry 
at the interface of B-C junctions, which could be shown 
to be independent of the periodic boundary conditions. 

Finally, we hope that further studies of suitable model 
systems similar to the structures presented in this study 
will shed some light on the formation and stability of 
other nanotubular heterojunctions and networks. Ac- 
cording to our studies nanotubular carbon and boron ma- 
terials seem to be compatible, and given the large number 
of structurally related materials (see introduction), our 
results also suggest the possibility of a large variety of 
novel heterogeneous nanotubular materials forming ex- 
tended network systems. By taking advantage of the 
fortunate structural and electronic properties of hetero- 
geneous nanotubular systems, one might even hope for 
the development of novel compound nanotubular devices 
in the future. 
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